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In this short letter we show that the excess of events in the decay of Higgs to two photons reported 
, by ATLAS and CMS can be easily accommodated in a flavor renormalizable three Higgs doublet 

model (3HDM). The model is consistent with all fermion masses, mixing angles, and flavor changing 

■ neutral current constraints. 

(N 

<^ , Recent results by ATLAS and CMS suggest an excess in the decay — >■ 77 with respect to the Standard Model 
D ' (SM) expectation [l|, 0]. Since the coupling of the Higgs to photons is radiative one possible reason for the excess, 
Mh , without changing the rest of the decays, is the existence of new charged particles with masses of hundreds of GeV. 

■ On the other hand there does not seem to be any significant deviation in the giuon fusion rate nor in the decay rates 
of the Higgs into gauge bosons, therefore the new particles responsible for the excess on the rate to photons have to 
be colorless and preferably not spin 1, in order to avoid any mixing with the W or the Z. Thus it seems that the 
existence of one or more spin or spin 1/2 charged particles could explain the diphoton rate measured at the LHC. 

I , In this short letter we show that it is possible to accommodate the excess in the diphoton rate within a renormalizable 
three Higgs doublet model (3HDM) whose motivation is to describe, in a minimal way, the flavor structure of the 
SM. Having extra Higgs doublets imply, upon electroweak symmetry breaking, the existence of new charged spin zero 
particles, two in this particular model. We are going to show that it is possible to satisfy all flavor constrains in this 
model as well as to generate an extra contribution to the decay rate of the Higgs to photons due to the two new 
charged degrees of freedom. 

The model: We explore a renormalizable flavor model based on the cyclic group Z5 that contains three Higgs 
. doublets. We chose Z5 as it has been shown to be the smallest Abelian symmetry that in the context of 3HDM 
' leads to the Nearest Neighbour Interaction (NNI) Yukawa textures in the quark sector Furthermore, the lepton 
T-H ■ sector of the model completely reproduces the one described in Q, where no right-handed neutrinos arc introduced 
' and neutrino Majorana masses are generated radiatively through the presence of an SU(2) singlet field charged under 
I both Hypercharge and Lepton number. 
^f) • The model particle content consists of the Standard Model fermion fields, three SU(2) doublets (Hypercharge 
y = 1/2) called H ~ [H, $1, $2), and a singlet scalar field 77 with Hypercharge y ~ —1 and Lepton number L — 2. 
The Z5 charge assignments for the fermion fields are given by 

^Ll ~ (0,-3,-1), Ufl ~dfl~efl~(3,0,2), (1) 

where Ql and Ll denote the quark and lepton left-handed SU(2) doublet SM fields respectively, while ur, dpi, and 
denote the quark and charged lepton right-handed SU(2) singlet SM fields. The scalar sector assignments are given 
by 

H = (if, $1, $2) ^(0,-1,1), ry~(-l). (2) 

These assignments lead to NNI textures for the quarks and charged lepton Yukawa matrices, while the neutrino 
mass matrix is of the form 
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The scalar potential for % is given by (see Eq. (12) of Q) 
+ A3 {<b\H^\H + /i.c) , 



(4) 



where a = 1,2 and the terms proportional to and ^12 are soft breaking terms required in order to obtain 
the correct electromagnetic invariant vacuum. The scalar 77 in the model plays its role in the neutrino sector where 
it provides the necessary Lepton number violation in order to generate Majorana masses. Due to its heavy mass 
and small mixing, it does not play a role in the phenomenology of the charged scalars present in Ti. When the 
diagonalization of the scalar sector is performed we thus only consider the potential in Eq. (|4]) . 

As discussed in 0, |3| the textures obtained in this model can reproduce all the observed masses and mixing angles 
in the quark and charged lepton number sectors, as well as the squared mass differences (for inverted hierarchy only) 
and mixing angles in the neutrino sector. Furthermore, it was also shown in that the model is able to satisfy the 
strong flavor changing neutral current constraints coming from K — K mixing. 

Results: A numerical scan of the parameter space of the model was performed in order to find regions where the 
model successfully reproduces/satisfies i) all quark masses and mixing angles, ii) charged lepton masses, squared 
neutrino mass differences and lepton mixing angles, iii) FCNC bounds from K — K mixing. 

The regions that satisfy all these constraints are then used to determine the contribution to the decay of the lightest 
neutral CP-even scalar, hP into two photons. This decay is computed using the expression Q 
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where i runs over scalars, fermions and vector bosons in the loop with charge e; and color factor Nc, and where the 
Fi factors are given by 



^^0 = [t(1 - r(r))] Gabc 

Fi/2 = [-2r(l + (l-T)/(T))]Ga„c 

Fi = [2 + 3r + 3r(2-T)/(r)]Gaf,c 



with T — (2mi/m/f)^ 



(6) 
(7) 
(8) 



sin-i(y(T77) 



if r > 1 



if r < 1 



(9) 



and where the factors Gate contain the model dependent information on mixing angles and couplings among fields a, 
b and c in the diagram vertices (Gabc = 1 V a, 5, c for the SM). 
Following 01 we define the ratio 



BR(/iO ^ 77) 



MODEL 



BR(/lO ^ 77) 



SM 



and take the best value for R~^^ obtained by combining in average the measurements presented in P, [2, |^, Q 



R^^ = 2.1 ±0.5 



(10) 



(11) 



The most recent results from ATLAS [10| quote a somewhat smaller ratio given by R^^ = 1.8 ± 0.3 with nih = 
126.6 GeV, whereas CMS has not presented any update in this channel. 

Figure [1] shows the values obtained for R^^ as a function of the vacuum expectation value oi H - vq - for three 
different sets of the parameters in the potential that satisfy the above conditions. Sets 1 and 2 correspond to cases 
where there are solutions that fall within a cr of the best value in (|Tl]) (shaded region with horizontal gray lines. The 
vertical green lines shade the region corresponding to the new ATLAS results) for some values of vq (Set 2 contains 
some points where the value of Rj^ falls outside). Set 3 on the other hand corresponds to a case where the model 
can give the same contribution than the SM and this happens for values of vq 233 — 236 GeV. The specific values 
for the parameter in the potential for each of the shown sets are given in Table H] Note that in particular Set 1 falls 
within the intersection of the two regions. 
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Ai 


Ao 


Aoi 


Ao2 


Al2 


Al2 


A21 


-^01 


-^02 


A3 A2 


Set 1 


-0.3782 


0.0721 


0.9473 


-0.2372 


-0.1784 


0.5587 


0.2107 


0.2954 


0.4455 


0.2098 0.3632 


Set 2 


0.9429 


-0.1628 


0.7062 


-0.3952 


-0.6729 


-0.0983 


-0.2106 


0.3305 


0.2143 


0.5857 0.4335 


Set 3 


-0.3169 


0.2017 


-0.9434 


0.4305 


0.9398 


0.2024 


-0.6206 


-0.1603 


0.8452 


-0.0159 0.7997 



Table I: Three sets of parameter values consistent with current experimental data on fcrmion masses and mixing angles and 
FCNC contributions to K — K mixing. The ratio R~,j obtained from these parameters as a function of the vev vq of H are shown 
in Figure [1] The vevs vi and «2 are obtained by requiring the hghtest neutral CP-even scalar to have a mass of = 125 GeV 
and from the relation vf — (246)^ as vq is varied in Figure[T] The soft braking parameters are taken to be /ifa = —(350)^, 
/iOl^ = -(400)^ and ^02'= -(450)^ 
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Figure 1: Contribution to i?-y^ as a function of vq for three different sets of parameters. All points in the plot correspond to 
cases that reproduce all fermion masses and mixing angles as well as FCNC bounds from K — K mixing. The specific values 
for the parameters in the potential are given for each set are given in Table m The shaded regions correspond to the 1 a range 
of R-^~i as given by (|lip (gray horizontal lines) and the latest ATLAS results (green vertical lines). 



In this letter we considered the recent diphoton rate excess presented by ATLAS and CMS in the context of a 
renormalizable flavor model with three Higgs doublets and based on Z5. Wc found that the model can accommodate 
the observed excess as well as all flavor observables including constraints from FCNC. 
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